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Abstract. Relativistic Faddeev equations for three-body scattering are solved 



^ I at arbitrary energies in terms of momentum vectors without employing a 



partial wave decomposition. Relativistic invariance is incorporated withing 
the framework of Poincare invariant quantum mechanics. Based on a Malfliet- 
Tjon interaction, observables for elastic and breakup scattering are calculated 



' and compared to non-relativistic ones. 
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A consistent treatment of intermediate energy reactions requires a Poincare 
symmetric quantum theory [T]. In addition, the standard partial wave decompo- 
sition, successfully applied below the pion-production threshold [2], is no longer 
OO . an adequate numerical scheme due to the proliferation of the number of par- 

tial waves. Thus, the intermediate energy regime is a new territory for few-body 
calculations, which waits to be explored. 
^ . This work addresses two aspects in this list of challenges: exact Poincare 

^ I invariance and calculations using vector variables instead of partial waves. In 

Ref. [3] the non-relativistic Faddeev equations were solved directly as function 
of vector variables for scattering at intermediate energies. A key advantage of 
this formulation lies in its applicability at higher energies, where the number of 
partial waves proliferates. The Faddeev equation, based on a Poincare invariant 
mass operator, has been formulated in detail in [3] and has both kinematical and 
dynamical differences with respect to the corresponding non-relativistic equation. 

The formulation of the theory is given in a representation of Poincare invari- 
ant quantum mechanics where the interactions are invariant with respect to kine- 
matic translations and rotations [5]. The model Hilbert space is a three- nucleon 
Hilbert space (thus not allowing for absorptive processes). The method used to 
introduce the NN interactions in the unitary representation of the Poincare group 
allows to input of e.g. high-precision NN interactions in a way that reproduces 
the measured two-body observables. However in this study we use a simpler in- 
teraction consisting of a superposition of an attractive and a repulsive Yukawa 
interaction with parameters chosen such that a bound state at Ed = -2.23 MeV 
is supported [IJ. Poincare invariance and S-matrix cluster properties dictate how 




Figure 1. The total elastic cm. cross section for elastic (top left) and for breakup scattering 
(bottom left) calculated from a Malfliet-Tjon type potential as function of the projectile kinetic 
laboratory energy. The labels 'R' ('NR') stand for relativistic (non-relativistic) calculations. The 
Faddeev calculations in the first order in t are marked with '1st', the converged full Faddeev 
calculations with 'full'. To show the difference, the percentage difference between the relativistic 
and corresponding non-relativistic calculations are displayed on the right. 



the two-body interactions must be embedded in the three-body dynamical gener- 
ators. Scattering observables are calculated using Faddeev equations formulated 
with the mass Casimir operator (rest Hamiltonian) constructed from these gen- 
erators. 

To obtain a valid estimate of the size of relativistic effects, it is important that 
the interactions employed in the relativistic and non-relativistic calculations are 
phase-shift equivalent. We follow here the suggestion by Coester, Piper, and Ser- 
duke (CPS) in constructing a phase equivalent interaction from a non-relativistic 
2N interaction [6] by adding the interaction to the square of the mass operator. 
In this CPS method the relativistic interaction can not be analytically calculated 
from the non-relativistic one. However, there is a simple analytic connection be- 
tween the relativistic and non-relativistic two-body t-matrices 

t.e(p,p';2i?-') = ^ == t„.(p,p';2ii;7), (1) 



where 2Ep^^ = 2\J m? + and 2Ep'^ ~ m ~'~ ^"^^ This relativistic two-body t- 
matrix tre(p, p'; 2Ep'^^) is scattering equivalent to the non-relativistic one at the 
same relative momentum p [7|. This t-matrix serves then as input to obtain the 
Poincare invariant transition amplitude of the 2N subsystem embedded in the 
three-particle Hilbert space via a first resolvent method as layed out in Ref . [1] . 
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Figure 2. The exclusive differential cross section for the reaction ^H(p,2p)n at 508 MeV labo- 
ratory projectile energy for different proton angle pairs 9i — 02 symmetric around the beam axis 
as function of the laboratory kinetic energy of one of the outgoing protons. The meaning of the 
curves are the same as in Fig.[TJ except that here '1' denotes the 1st order Faddeev calculation, 
'F' the fully converged one. In the curves labeled -Rkin only relativistic kinematics is taken into 
account. The data are taken from Ref. |10] . 



By construction, differences in the relativistic and non-relativistic calcula- 
tions first appear in the three-body calculations. Those differences are in the 
choice of kinematic variables (Jacobi momenta are constructed using Lorentz 
boosts rather then Galilean boosts) and in the embedding of the two-body in- 
teractions in the three-body problem, which is a consequence of the non-linear 
relation between the two and three-body mass operators. These differences mod- 
ify the permutation operators and the off-shell properties of the kernel of the 
Faddeev equations [9]. 

In Fig. [1] the total cross sections for elastic and breakup cross sections are 
displayed as function of the projectile kinetic energy up to 1.5 GeV obtained from 
our fully converged relativistic Faddeev calculation as well as the one obtained 
from the first-order term, T^** = tP, with P being the permutation operator for 
three identical particles. It is obvious that, especially for energies below 300 MeV, 
the contribution of rescattering terms is huge. However, for extracting the size of 
relativistic effect, it is more useful to consider the relative difference between the 
relativistic and non-relativistic calculations. In first order, there is essentially no 
effect in the total elastic cross section, which is consistent with the observation 
that when the relativistic two-body i-matrix is constructed to be phase-shift 
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equivalent to the non-relativistic one. Making the same comparison with fully 
converged Faddeev calculations indicates that relativistic effects in the three- 
body problem increase the total cross section for elastic scattering with increasing 
energy, whereas it is slightly reduced in the total breakup cross section. 

Considering exclusive breakup reactions, differences between a relativistic 
and non-relativistic calculation can be more pronounced and strongly depend 
on the specific configuration. Though our two-body force is simple, we compare 
to a ^H(p,2p)n experiment at 508 MeV [10] to see if our calculation captures 
essential features of the measurement. Differences in the predictions of our rel- 
ativistic and non-relativistic calculations are very pronounced at this energy as 
can be seen in Fig. [21 which shows selected angle pairs 9i — O2 from Ref. [lOj . 
which are symmetric around the beam axis. The cross section is plotted against 
the laboratory kinetic energy of one of the outgoing protons. It is interesting to 
observe that for smaller angle pairs the relativistic cross sections (RF) are con- 
siderably larger than the non-relativistic ones (NRF). For larger angle pairs the 
situation reverses. It is further noteworthy, that in the configurations of Fig. [21 
which are close to quasi-free, rescattering effects (or equivalently higher order 
contributions of the Faddeev multiple scattering series) are very small. To show 
that peak-positions are given by kinematics, we added a curve labeled 'Rfcm' to 
the figures, which stands for a non-relativistic calculation in which only kine- 
matics and phase space factors are replaced by the relativistic ones. We want to 
note that the above comparisons do not involve a non-relativistic limit, instead 
relativistic and non-relativistic three-body calculations with interactions that are 
fit to the same two-body data are compared. All of the differences are due to the 
different ways two-body dynamics is incorporated in the three-body problem. 
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